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Memory is a hallmark of adaptive immunity, wherein lymphocytes mount a superior response to a previously en- 
countered antigen. It has been speculated that epigenetic alterations in memory lymphocytes contribute to their func- 
tional distinction from their naive counterparts. However, the nature and extent of epigenetic alterations in memory 
compartments remain poorly characterized. Here we profile the DNA methylome and the transcriptome of B-Iymphocyte 
subsets representing stages of the humoral immune response before and after antigen exposure in vivo from multiple 
humans. A significant percentage of activation-induced losses of DNA methylation mapped to transcription factor 
binding sites. An additional class of demethylated loci mapped to Alu elements across the genome and accompanied 
repression of DNA methyltransferase 3A. The activation-dependent DNA methylation changes were largely retained in 
the progeny of activated B cells, generating a similar epigenetic signature in downstream memory B cells and plasma cells 
with distinct transcriptional programs. These findings provide insights into the methylation dynamics of the genome 
during cellular differentiation in an immune response. 



[Supplemental material is available for this article.] 

Epigenetic modifications play important roles in regulating cellu- 
lar differentiation events. One such epigenetic modification, DNA 
methylation, occurs on cytosine residues primarily at CpG di- 
nucleotides in mammals. The role of DNA methylation in regu- 
lating cellular differentiation from pluripotent and multipotent 
progenitors has been demonstrated through functional analysis of 
animals deficient in DNA methyltransferases (DNMTs) (Li et al. 
1992; Okano et al. 1999; Tadokoro et al. 2007; Broske et al. 2009; 
Wu et al. 2010), as well as from recent genome-wide studies com- 
paring the DNA methylome of various differentiated cell types and 
their precursors (Meissner et al. 2008; Lister et al. 2009; Ji et al. 
2010; Hodges et al. 2011; Bock et al. 2012). In the context of the 
immune system, mutations in the DNMT3B gene are causal for 
the development of ICF syndrome (immunodeficiency, centro- 
mere instability, and facial anomalies syndrome), a rare autosomal 
recessive immune disorder (Hansen et al. 1999; Xu et al. 1999). 
Despite having a normal number of mature B cells, ICF patients 
lack memory B cells as well as plasma cells (PCs) (Blanco-Betancourt 
et al. 2004), suggesting the involvement of DNMT3B and possi- 
bly of DNA methylation in regulating late stages of lymphocyte 
maturation. 

Upon activation by antigenic stimulation in a T-cell-dependent 
B-cell immune response, naive B cells enter the germinal-center 
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(GC) reaction in secondary lymphoid organs. Within GCs, B cells 
activated by antigenic stimuli clonally expand and their immu- 
noglobulin (Ig) gene loci are subjected to somatic hypermutation 
and class-switch recombination (Victora and Nussenzweig 2012). 
These genetic alterations are critical for the maturation of GC B 
cells to post-GC cell types that subsequently produce high-affinity 
antibodies against foreign pathogens. Upon exiting the GC, B cells 
either differentiate into antibody-producing long-lived PCs or 
alternatively become memory B cells that provide long-term im- 
munity against the same pathogen (Shapiro-Shelef and Calame 
2005). During a secondary immune challenge, memory B cells 
more rapidly undergo a proliferative burst, and then differen- 
tiate into PCs in a facilitated manner compared to naive B cells 
(McHeyzer- Williams and McHeyzer- Williams 2005; Lanzavecchia 
and Sallusto 2009). Compared to naive B cells, the memory coun- 
terparts express B-cell receptors with higher affinity to the same 
antigen (Pascual et al. 1994), constitutively express costimulatory 
molecules on their cell surface (Liu et al. 1995), and have lower 
expression of transcription factors (TFs) important for maintaining 
cellular quiescence (Good and Tangye 2007). These unique fea- 
tures decrease the threshold of activation in memory B cells and 
allow them to quickly enter the cell cycle upon restimulation. 
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Aside from these key differences, naive and memory B cells possess 
highly similar gene expression programs (Klein et al. 2003), and it 
remains unclear how memory B cells can more efficiently repro- 
gram their transcriptional profiles to specify a PC fate. It has been 
speculated that epigenetic alterations in naive and memory lym- 
phocytes contribute to their functional outcomes (Messi et al. 
2003; Kersh et al. 2006; Cuddapah et al. 2010). Nonetheless, the 
degree of epigenetic differences in these two cell types remains 
undefined. It is also unclear whether DNA methylation plays a role 
in specifying an effector vs. a memory cell fate in lymphocytes 
during a humoral immune response. 

The global methylation landscape of the total B-cell fraction 
in peripheral blood was previously characterized, revealing distri- 
bution of this epigenetic mark at different genomic features (Rauch 
et al. 2009). To further understand the dynamics of DNA methyl- 
ation changes during a B-cell immune response, naive, GC, mem- 
ory and PC populations were purified ex vivo from inflamed tonsils 
of eight individual humans for global DNA methylation analysis. 
We observed widespread DNA methylation changes between naive 
and GC B cells. In contrast, fewer changes in DNA methylation 
occur upon resolution of activated, GC B cells into either memory 
or PC alternative cell fates. These results demonstrate that memory 
B cells and PCs share a highly similar methylome, despite having 
very distinct transcriptional programs and that naive and memory 
B cells, which have highly similar transcriptomes, have vastly 
different methylomes. In agreement with genome-scale studies in 
other cell types (Stadler et al. 201 1; Bock et al. 2012; Thurman et al. 
2012), differentially methylated regions (DMRs) in GC B cells were 
enriched for TF binding sites or regulatory elements. A second class 
of loss-of -methylation events localized to Alu elements, consistent 
with the preferential targeting of DNMT3 family de novo meth- 
yltransferases to repetitive elements in the genome (Jeong et al. 
2009). The methylation dynamics dXAlu elements suggests a reg- 
ulatory role of the repetitive fraction of the genome in the GC 
transcriptional program and/or its associated chromatin architec- 
ture, and a potential contribution to genomic instability in GC B 
cells (Ranuncolo et al. 2007). We propose that global reprogram- 
ming of the epigenome during the GC reaction represents an event 
that permits an adult somatic cell to differentiate into multiple, 
diverse cell types. In addition, the epigenetic signature of memory 
B cells may contribute to their superior ability, compared to that 
of naive B cells, to rapidly differentiate into PCs upon antigen 
rechallenge during a secondary immune response. 

Results 

Widespread alteration of the DNA methylation landscape 
in activated B cells is conserved in post-GC cell types 

Naive, GC, memory B cells, and PCs were purified from inflamed 
tonsils by FACS (Fig. 1A,B). The four cell types can be distinguished 
by IgD and CD38 cell surface expression within CD20 + B cells 
(Pascual et al. 1994; Klein et al. 2003). In monitoring expression of 
key molecules by quantitative RT-PCR to validate cell purity we 
observed significant changes in expression of all three mamma- 
lian DNMTs (DNMT1 , DNMT3A, and DNMT3B) (Jones and Liang 
2009) in B-cell subsets during the course of an immune response 
(Fig. 1C). While the expression level of DNMT3A dramatically 
decreased upon activation by antigen (naive to GC B-cell transi- 
tion), DNMT1 and DNMT3B were up-regulated at this stage, 
consistent with the role of DNMT1 in replication-dependent 
maintenance methylation (Leonhardt et al. 1992) and with the 



role of DNMT3B in late-stage B-cell differentiation (Blanco- 
Betancourt et al. 2004). Post-GC memory B cells and PCs rever- 
ted the expression pattern of DNMTs to that of their naive 
precursors (Fig. 1C). Immunohistochemical staining of tonsil sec- 
tions confirmed robust expression of DNMT1 and minimal ex- 
pression of DNMT3A in GC B cells (Fig. ID). In contrast, regions 
surrounding the GC comprised primarily of naive B cells have the 
opposing expression pattern of the corresponding DNMT. Up- 
regulation of DNMT3B in GC B cells at the protein level was also 
confirmed by immunofluorescence (Fig. IE). Expression changes in 
all three forms of DNMTs upon immune activation pointed 
to possible reprogramming of the DNA methylation landscape in 
GCB cells. 

To examine the dynamics of DNA methylation during the 
immune response, global DNA methylation profiling was per- 
formed on the four tonsillar B-cell subsets from eight individuals 
using the methylated CpG island recovery assay (MIRA) (Rauch 
et al. 2006) and analyzed on microarrays (MIRA-chip). While the 
number and genomic coverage of autosomal regions enriched in 
DNA methylation peaks varied by cell type (Supplemental Table SI) 
and also by subject (data not shown), DNA methylation dynamics 
relevant to immune activation and the subsequent cell differenti- 
ation programs were presumed to be consistent across all indi- 
viduals. DMRs were identified for each pair of cell types across all 
subjects to assess the DNA methylation changes during the im- 
mune response (Supplemental Table S2). There were a greater 
number of DMRs associated with immune activation, at the tran- 
sition from naive to GC B-cell stage (Fig. 2A; Supplemental Table S2), 
than with any other developmental transition. Notably, immune 
activation-induced DMRs were dominated by loss-of-methylation 
events (Fig. 2A), in agreement with a previous study profiling the 
status of 50,000 CpGs during activation (Shaknovich et al. 2011). 
Substantially fewer DMRs were identified between GC B cells 
and memory B cells or PCs (Fig. 2A; Supplemental Table S2). DNA 
methylation levels at activation-induced DMRs were highly similar 
when comparing GC, memory B cells, and PCs (Fig. 2B, left panel). 
Principle component analysis (PCA) also indicated that the 
methylation signature of GC and post-GC cell types cluster closely 
together, clearly distinct from naive B cells (Fig. 2B, right panel). 
These data indicate that the majority of activation-induced DNA 
methylation changes in the GC are inherited by both memory 
B cells and PCs. Loss-of-methylation events in GC B cells were 
equally distributed near promoter-proximal regions and promoter- 
distal regions (Fig. 2C,D). In contrast, gain-of-methylation events 
were overrepresented in promoter-proximal regions (Fig. 2C,D). 
The genomic distribution between the two types of DMRs is sig- 
nificantly different by the \ 2 test, P < 0.0001, indicating the pos- 
sibility of distinct regulatory roles of gain- and loss-of-methylation 
events. 

Widespread loss-of-methylation could be a result of passive or 
active demethylation (Wu and Zhang 2010). Two major active 
demethylation pathways have been established in mammals: (1) 
activation-induced cytosine deaminase (AICDA)/APOBEC-dependent 
deamination of 5-methyl-cytosine (5mC) followed by base ex- 
cision repair; and (2) oxidative demethylation by TET proteins. 
Although AICDA is highly expressed in the GC, recent genome- 
wide analysis of AICDA occupancy indicates a strong pref- 
erence for promoter-proximal regions (Yamane et al. 2011). 
However, we did not observe a biased distribution of loss-of- 
methylation DMRs near TSSs in GC B cells (Fig. 2C,D). Further, 
AICDA/APOBEC enzymes have substantially reduced activity 
on 5mC relative to unmodified cytosine (Nabel et al. 2012). 
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Figure 1. Expression pattern of DNMTs in tonsillar B-cell subsets. (A) Cartoon illustrating the stages of B-cell development and T-dependent B-cell 
immune response. B cells are differentiated from hematopoietic stem cells (HSCs) in the bone marrow and released into the periphery upon maturation. 
Naive B cells encounter antigen-primed T cells in secondary lymphoid organs and enter the germinal center (GC) reaction followed by differentiation into 
memory B cells or antibody-secreting plasma cells (PCs). (B) Purification of B-cell subsets. CD20, IgD, and CD38 cell surface expression were used to distinguish 
the cell populations: Naive B cells (CD20 + lgD + CD38'°), GC B cells (CD20 + lgDCD38 + ), memory B cells (CD20 + lgDCD38~), and PCs (CD20 lo/+ lgD~CD38 hi ). 
The percentage of each subset within tonsillar B cells is shown in the FACS plot in the middle panel. Naive, GC, memory, and PC populations were purified 
from eight individuals for gene expression and DNA methylation profiling on microarrays. The purity of each cell population (>90%) was confirmed by 
flow cytometry analysis of post-sorted cells (right panel). (C) Transcript abundance of DNMTs in B-cell subsets by quantitative RT-PCR. DNMT expression 
levels are normalized to actin expression level in each cell type. The expression level in one biological replicate of naive B cells is arbitrarily set to one. Error 
bars denote standard deviation of three biological replicates. (D) Expression patterns of DNMT1 and DNMT3A in tonsils. (Upper left, upper right) Im- 
munohistochemical staining of DNMT1 at 200x and 400x magnifications, respectively. (Lower left, lower right) Immunohistochemical staining of 
DNMT3A at 200x and 400x magnifications, respectively. (£) Up-regulation of DNMT3B in GC B cells. Tonsil sections were costained with anti-DNMT3B 
(red) and anti-MTA3 (green). Nuclei were labeled with DAPI (blue). DNMT3B expression colocalized with MTA3, which is a marker for GC B cells (Fujita 
et al. 2004). Scale bar, 50 |xm. 



Therefore, AICDA may not be the mediator of the alterations 
observed here. 

To explore the TET pathway in this system, we examined the 
global level of 5hmC, the product of TET proteins, in naive and GC 
B cells by immunofluorescence. While the level of 5mC remained 
unchanged in the two B-cell subsets (Fig. 2E), we observed a de- 



crease in 5hmC in GC B cells (Fig. 2F), arguing against the role of 
this demethylation pathway in the GC. The most likely explana- 
tion for the site-specific DNA methylation changes observed in GCs 
is passive demethylation rather than an active process. We envision 
two possible explanations that are not mutually exclusive: (1) De- 
clines in DNMT3A levels coupled with a massive burst of pro- 
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Figure 2. Widespread alteration of DNA methylation in B-cell subsets. (A) Assessment of global DNA methylation changes during immune activation 
(naive [N] vs. GC) and differentiation of GC B cells into memory or PCs (GC vs. memory [Mem] and GC vs. PC) by MIRA-chip analysis. The bar graph 
displays the number of DMRs with increased methylation (gain of 5mC) in red or decreased methylation (loss of 5mC) in green in each stage transition. (B) 
The methylation score is computed for each DMR between naive and GC B cells in all MIRA-chip samples and displayed in a heat map (left panel). Each row 
represents a DMR between naive and GC. Each column represents an individual MIRA-chip sample. All eight samples from each cell type clustered together 
in the heat map. Cell types are represented by different colors, as indicated in diagram beside the heat map. PCA analysis of methylation score of DMRs for 
all cell types is shown in the right panel. Each circle represents an individual biological sample, x- and /-axes represent first and second principal com- 
ponents, respectively. The degree of variance in each principal component is indicated in parentheses. (C) Frequency and genomic distribution of DMRs 
between naive and GC B cells. The bar graphs display the number of DMRs that are <1 0 kb or >1 0 kb from the nearest annotated TSS. The left panel 
represents DMRs with loss of methylation between naive and GC B cells (N > GC), while the right panel represents DMRs with gain of methylation between 
naive and GC B cells (N < GC). (D) Distribution of DMRs relative to the nearest TSS. The frequency of both types of DMRs, loss-of-methylation DMRs (N > 
GC), and gain-of-methylation DMRs (N < GC), in 1 -kb increments from the TSS is shown in the graph. (£) Distribution of methylcytosine (5mC) in tonsils. A 
tonsil section was double stained for MTA3 to demarcate GC B lymphocytes (green) and 5mC (red). Nuclei were labeled with DAPI (blue). There was no 
detectable difference in the global level of 5mC in MTA3 + GC B cells and in non-GC B cells. Scale bar, 50 |jim. (F) Distribution of hydroxy-methylcytosine 
(5hmC) in tonsils. A tonsil section was double stained for MTA3 to demarcate GC B lymphocytes (green) and 5hmC (red). Nuclei were labeled with DAPI 
(blue). The global level of 5hmC decreased in in MTA3 + GC B cells compared to non-GC B cells. Scale bar, 50 ixm. 
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liferation lead to a failure to maintain DNA methylation levels; and 
(2) global alteration in nuclear architecture and nucleosome posi- 
tioning results in a change in substrate quality for DNMTs. 

B-Iymphocyte activation establishes DNA methylation 
reprogramming that poises memory B cells for recall responses 

To gain insight into the relationship between DNA methylation 
and gene expression, we performed global transcriptional profiling 
in the purified cell populations used for DNA methylation analysis. 



Selected genes diagnostic of the developmental transitions occur- 
ring during the immune response were manually validated by 
RT-PCR and inspected from the microarray data (Supplemental 
Fig. S1A ; B). Unsupervised hierarchical clustering and principal 
component analysis (PCA) of the gene expression profiles in- 
dicated that naive and memory B cells share a highly similar 
transcriptional program, but are distinct from GC as well as PC cell 
types (Fig. 3, far left panels). More than 1200 genes were differ- 
entially expressed between each pair of cell types examined, the 
exception being the comparison between naive and memory cells 
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Figure 3. B- lymphocyte activation establishes DNA methylation reprogramming that poises memory B cells for recall responses. (Top) Heat maps 
displaying hierarchical clustering of genes differentially expressed (ANOVA P < 0.01 and fold change > 3) and genomic features with differential meth- 
ylation scores (ANOVA P < 0.01). The total number of genes or genomic features is indicated. (Middle) PCA analysis of the corresponding differentially 
expressed genes or differentially methylated genomic features. Each circle represents an individual MIRA-chip sample; cell types are represented by 
different colors, x- and /-axes represent first and second principal components, respectively. The degree of variance in each principal component is 
indicated in parentheses. (Bottom) Analysis of genomic features by average score in naive B cells (N) vs. PC; naive B cells (N) vs. memory B cells (M); and 
memory B cells (M) vs. PC. The bar graph displays the number of genomic features with increased methylation (up) in red and decreased methylation 
(down) in green. 
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where only 158 genes were differentially expressed (Supplemental 
Fig. SIC). 

The extent of DNA methylation difference between the four 
B-cell subsets was further investigated using a mixed-effect model 
to analyze the differences in average methylation score at various 
defined genomic features (Weber et al. 2007). Strikingly the close 
similarity of memory and naive cell populations in terms of tran- 
scriptional profile contrasts sharply with the pattern of DNA meth- 
ylation at all genomic features examined. Although there were dif- 
ferentially methylated genomic features between memory B cells 
and PCs, these two cell types clustered closely at most features 
(Alus, CpG islands, and CpG island shores) (Fig. 3, top panel), sug- 
gesting a high degree of similarity between the two methylomes. 
The first principal component of PCA resolved GC B, memory B, 
and PC from naive B cells for all four genomic features, while a 
much smaller degree of variance was detected between GC, memory 
B, and PCs (Fig. 3, middle panel). The genomic feature analysis was 
consistent with the DMR analysis (Fig. 2B; Supplemental Table S2), 
suggesting that DNA methylation is globally altered upon activa- 
tion in naive B cells, while the dynamics are more subtle as cells 
partition between memory B-cell vs. PC compartments as their 
post-GC fate. In agreement with hierarchical clustering and PCA, 
there are significantly fewer differentially methylated genomic 
features between memory B cells and PCs, compared to naive vs. 
PCs or naive vs. memory B cells (Fig. 3, bottom panel; Supple- 
mental Tables S2, S3). These data demonstrate that antigen-driven 
activation in B cells induces global DNA methylation reprogram- 
ming, resulting in a common epigenetic signature in memory B 
cells and PCs, despite the adoption of very different transcriptional 
programs in post-GC cell types. 

To determine whether the DNA methylation profiles described 
here may contribute to the functional superiority of memory cells 
in differentiation into PCs, we investigated whether genomic re- 
gions differentially methylated between naive and memory B cells 
are associated with genes or regulatory elements with known 
functions in immune responses. DMRs between naive and mem- 
ory B cells that were associated with genes (within 10 kb of TSS) were 
divided into two categories: naive > memory or naive < memory in 
methylation signal. Both categories were enriched for genes with 
functions associated with immune response and hematological 
functions, gene expression, cell growth and proliferation, and cel- 
lular functions integral to lymphocyte activation (Supplemental 
Table S4). Genomic feature analysis at promoters (-1 to +1 kb from 
TSS) and CpG island shores (located within 10 kb of gene pro- 
moters) revealed highly similar enrichment in functions relevant to 
immune response (data not shown). Exemplar genes having im- 
mune-related functions that differ in methylation pattern (Supple- 
mental Fig. S2A) include critical regulators of B-cell development 
(£4X5), of immune responses (RELA), and of cell death/survival 
(BCL2L1 and RUNX3). Importantly, the expression level of these 
genes, while similar in naive and memory B cells, changes signifi- 
cantly upon differentiation into PCs (Supplemental Fig. S2B). The 
association of DNA methylation changes with elements associated 
with immune activation and the immune response supports the con- 
cept that the memory B-cell epigenome is poised to facilitate a more 
rapid and robust activation response than their naive counterparts. 

Enrichment of regulatory elements and TF binding sites 
at activation-induced DMRs 

To further investigate the function of DNA methylation changes 
associated with immune activation, we first determined whether 



activation-induced DMRs were enriched for regulatory elements. 
The degree of overlap between DMRs and DNase I hypersensi- 
tive sites (DHSs) from ENCODE DNase I-seq data sets (Thurman 
et al. 2012) was queried. Both loss-of-methylation and gain-of- 
methylation DMRs were highly enriched for DNase I clusters, 
which consist of combined DHS identified across 125 diverse cell 
and tissue types (Fig. 4 A; Thurman et al. 2012). The majority of the 
DHS overlapping DMRs were present in primary B cells (Fig. 4A), 
consistent with a relationship between DNA methylation and 
regulatory factor occupancy patterns (Stadler et al. 201 1; Bock et al. 
2012; Thurman et al. 2012). We further searched for the presence of 
DNA binding motifs of several TFs important for B-cell development 
within the DMRs (Table 1). In line with the DNase I HS analysis, we 
identified several TFs, including MYC and BCL11A, with motifs 
that were significantly enriched within loss-of-methylation DMRs 
in GC B cells (Fig. 4B; Table 1). Analysis of MYC by ChlP-PCR 
(Supplemental Fig. S3A) and of MYC and BCL11A ChlP-seq data 
from ENCODE (The ENCODE Project Consortium 2011) in a 
B-lymphoblastoid cell line indicates that the DMRs frequently 
overlap with the binding sites of these TFs (Fig. 4B; Supplemental 
Fig. S3 A). Further strengthening the connection between DMRs 
and regulatory DNA, we observed increased levels of histone 3 
lysine 27 acetylation (H3K27ac), consistent with these regions 
being active enhancers (Creyghton et al. 2010), in GC B cells at 
several loss-of-methylation DMRs (Supplemental Fig. S3B). 

The relationship between DNA methylation dynamics and 
gene expression was further examined. Gene set enrichment 
analysis (GSEA) (Subramanian et al. 2005) revealed that genes with 
down-regulated transcription during the naive to GC or the naive 
to PC differentiation events were significantly enriched for those 
that have overall gained methylation in GC B cells (Figs. 4C; Sup- 
plemental Table S5). PAX5, a TF critical for B-lineage commitment 
whose down-regulation is required for the activation of the PC 
transcriptional program (Delogu et al. 2006), exemplifies genes 
gaining methylation on activation (Supplemental Fig. S2B). These 
results are consistent with the observation that gain-of-methyla- 
tion DMRs are preferentially enriched near transcription units and 
within genie regions (Fig. 2D). We did not observe a significant 
association between loss of DNA methylation on activation and 
gene expression, perhaps reflecting the likelihood that many of 
these methylation losses occur at distal regulatory DNA that is 
difficult to assign to a particular gene. In any case, these analyses 
indicate that DNA methylation dynamics in B-cell subsets repre- 
senting various stages of the immune response provide insights 
into the regulatory elements that govern the GC reaction and 
transcriptional program. 

Widespread loss of DNA methylation at Alu elements 
on activation 

About half of activation-induced DMRs overlapped with DHS 
(Fig. 4A), the mechanistic underpinning of the remaining DMRs 
was unclear. When mapping DMRs to various genomic features, 
a large proportion of DMRs between naive and GC B cells were 
adjacent to Alu elements (Fig. 5 A; Supplemental Table S2). This 
observation is contrary to the widely held view that methylation at 
repetitive elements is constitutive in normal somatic cell types 
(Hellmann-Blumberg et al. 1993). Loss-of-methylation DMRs are 
especially predominant over gain-of-methylation events near this 
genomic feature (Fig. 5 A). 

To explore the prevalence of loss-of-methylation at Alu ele- 
ments in GC B cells, we analyzed the average methylation score at 
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Figure 4. Activation-induced DMRs reveal c/'s- regulatory elements in GC B cells. (A) Overlap of DMRs 
with DNase I hypersensitive sites (DHSs). The percentage of loss-of-methylation (N > GC) or gain-of- 
methylation (N < GC) DMRs intersecting DNase I clusters from any cell type or DHS from primary CD20 + 
B cells is shown. (B) Example of a DMR overlapping a DNase I HS as well as MYC and BCL1 1 A binding 
sites in a lymphoblastoid cell line. The bar graphs display the average MIRA-chip signals from eight 
biological replicates in naive, GC, memory, and PC cells at the genomic region on chromosome 19, as 
indicated on the x-axis. The location of an annotated gene (TNFSF9) is shown below the bar graphs 
relative to the genomic coordinates on chromosome 1 9, and the orientation of the arrow indicates the 
direction of transcription. Black rectangles below the plots indicate the location of MYC and BCL1 1 A 
binding sites determined by ChlP-seq data from the GM1 2878 lymphoblastoid cell line (performed by 
Iyer's group at University of Texas-Austin and Myer's group at HudsonAlpha, respectively) as de- 
termined by the ENCODE project (The ENCODE Project Consortium 201 1). (C) Analysis of the asso- 
ciation between DNA methylation and gene expression changes by GSEA. The -log of enrichment 
P-value for each gene set (differentially expressed genes in each pair of cell types) is shown. Gene sets 
with -log of P-value above the dotted line (>2 or P< 0.01) are considered to be statistically significant. 



each Alu element covered in our experiments (Fig. 5A; Supple- 
mental Fig. S4). Since the microarrays lack coverage within re- 
petitive elements, we used probes adjacent to repeats (within 
250 bp) as a surrogate for their methylation status in subsequent 
analysis as these closely correlated with the methylation levels 
within candidate Alu elements assessed by genomic bisulfite 
sequencing (Supplemental Fig. S4). In agreement with the DMR 
analysis, differentially methylated Alus lost methylation in GC B 
cells and remain hypomethylated in post-GC cell types (Fig. 3). 
PCA analysis of differentially methylated Alus again indicated 
dissimilarity between naive B cells and all other cell types arising 
after immune activation (Fig. 3). 

We further investigated the dynamics of Alu methylation in 
naive and GC B cells at the base pair level, by performing next 
generation sequencing of bisulfite-converted MIRA-enriched DNA 
fragments from naive and GC B cells (MIRA-seq) (Serre et al. 2010). 
The AluY subfamily of elements had the highest sequencing cov- 
erage (Supplemental Table S6), reflective of the high CpG content 
within this Alu family (Supplemental Fig. S5A). Since the majority 
of MIRA-enriched CpG dinucleotides were methylated (Supple- 
mental Table S7), we identified differentially methylated Alus 



solely based on sequencing depth. Cross- 
analysis of MIRA-chip and sequencing 
data indicated strong correlation between 
the two data sets (Supplemental Fig. S5B; 
Supplemental Table S8). We observed 
similar numbers of gain- and loss-of- 
methylation events on activation in the 
older J and S Alu subfamilies, but loss-of- 
methylation events dominated at AluYs 
(Fig. 5B), the youngest class of Alu ele- 
ments. Roughly 10% of elements within 
the AluY family had methylation changes 
between naive and GC B cells, totaling over 
14,000 elements (Supplemental Table S6). 
Validation of differential methylation sta- 
tus in naive and GC B cells was performed 
in a panel of Alu elements by genomic bi- 
sulfite sequencing (Fig. 5C; Supplemental 
Table S9). The distinct methylation dy- 
namics within the AluY family suggest that 
either these events play a functional role 
in B-cell activation or that they reflect 
a global change in the nucleus occurring at 
this stage transition. Cross-analysis with 
ENCODE DNase I-seq data indicated that 
a substantially lower percentage of dif- 
ferentially methylated Alus, compared to 
DMRs, overlapped regulatory factor bind- 
ing sites (Supplemental Fig. S3C). While 
the genomic targets of retinoic acid re- 
ceptors (RARs) are frequently found in Alu 
elements (Laperriere et al. 2007), the pres- 
ence of the RAR motif was not enriched in 
differentially methylated A lus in GC B cells 
(Supplemental Fig. S3D), further suggesting 
that differential methylation of these ele- 
ments plays a biological role distinct from 
modulation of TF binding. 

Similar to all DMRs between naive 
and GC B cells, gain-of -methylation AluYs 
were found to be enriched near promoter- 
proximal regions and in transcription units (Fig. 6A). In contrast, 
loss-of-methylation AluYs were spread throughout the genome, 
mirroring the distribution of all elements within the AluY sub- 
family (Fig. 6A). Chromosome distribution analysis indicated loss- 
of-methylation events at AluYs were global and pervasive across 
the genome in GC B cells upon immune activation (Fig. 6B). 
Nonetheless, pathway analysis of genes associated with loss-of- 
methylation Alus indicated enrichment of genes with functional 
roles in cell growth and proliferation (Supplemental Table S10), 
a critical component of the GC reaction. Interestingly, while gain- 
of-methylation distributes evenly across the entire consensus 
Alu element in both Y and S subfamilies, loss-of-methylation 
events preferentially occurred at the 5' and 3' ends of the elements 
(Figs. 6C; Supplemental Fig. S5C,D). Since these regions have well- 
positioned nucleosomes (Englander et al. 1993; Englander and 
Howard 1995), the observed 5' and 3' loss-of-methylation regions 
could reflect repositioning of nucleosomes dXAlu elements in naive 
to GC B-cell transition. The widespread nature of Alu demethyla- 
tion may therefore be indicative of a global restructuring of chro- 
mosomal architecture and nucleosome positioning in B cells upon 
immune activation. 
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Table 1. Transcription factor motif analysis of DMRs between naive and GC B cells 

Naive > GC DMRs Naive < GC DMRs 



Transcription factor 


Motif sequence 


Enrichment P-value a 


Enrichment P-value a 


MYC 


CACGTG 


0 


0.1 1 


BCL11A(motif 1) 


CC[C|T][A|G]C 

L | J L | J 


0 


1.0 


RBPJ (RBP JK) 


CGTGGGAA 


8.22 X 10~ 7 


0.57 


EBF 


CCCNNGGG 


1.58 X 10" 7 


1.0 


TCF3 (E2A) 


CANNTG 


1.91 X 10 -5 


1.0 


FOS 


TGACTCA 


5.97 X 10 -4 


0.65 


SPI1 (PU.1) 


GGA[A|T] 


6.0 X 10 B 


1.53 x 10" 12 


NFKB1 


GGGACTTTCC 


0.08 


NA 


BCL11A(motif 2) 


GGCCGGAGG 


0.17 


NA 


STAT3 


TT[A|C]CA[A|G]GAA 


0.17 


0.82 


BCL6 


TC[C|T][A|T][A|C]GA 


0.25 


0.01 


POU2F1 (OCT1) 


ATTTGCAT 


0.34 


0.79 


STAT5 


TTCTCAGAA 


0.44 


0.51 


CTCF 


CC[A|T|G][C|G]NAG[G|A][G|T]GG[C|A|T][G|A|T][C|G] 


0.50 


0.47 


IKAROS 


[T|C]GGGAAA[A|T] 


0.95 


0.01 


PAX5 


GCAT[C|G]A[T|A]GCGT[A|G][C|A] 


NA 


NA 


SOX4 


ATACAAATGG 


NA 


NA 


IRF4 


AAAGGAAGTGAACCA 


NA 


NA 


PRDM1 


AAGTGAAAGTGAA 


NA 


NA 



The presence of B-cell-related TF motifs (Matthias and Rolink 2005) was searched within genomic regions covered by naive > GC or naive < GC DMRs. An 
enrichment P-value was obtained for each motif by comparing the observed number of motifs in the DMRs to the expected number of motifs in genomic 
regions (probes) covered by the microarray using the Poisson distribution. 
Significant P-values are defined as P < 0.01 (in bold). 



The genomic distribution and pervasiveness of differentially 
methylated AluYs prompted us to investigate whether the meth- 
ylation dynamics at Alu elements reflect a global change in DNA 
methylation. Bisulfite sequencing was performed in naive and GC 
B cells at a limited number of regions proximal to differentially 
methylated Alus, indicating that these neighboring regions were 
also undergoing methylation changes (Fig. 6D; Supplemental 
Table Sll). Since MIRA enrichment preferentially captures CpG- 
rich genomic regions (Supplemental Table S7), future whole- 
genome bisulfite sequencing in these cell types will more precisely 
reveal the extent of methylation change outside of Alu elements 
and at elements with lower CpG density. 

Analysis of DNA methylation content in genomic regions 
varying in gene density revealed an overall decreased methylation 
level upon B-cell activation in genomic regions with low gene 
density, while regions with high gene density conversely have 
increased methylation content (Fig. 6E). Since we did not observe 
a decrease in total 5mC content, this result indicates that there is 
a global redistribution of CpG methylation events at gene-rich vs. 
gene-poor regions. These data provide further support for a global 
restructuring of chromosomal architecture in B cells upon acti- 
vation. Collectively, the DNA methylation dynamics described 
here provide a useful platform to gain insights into the tran- 
scriptional network in a GC reaction, as well as uncovering po- 
tential novel regulatory elements in the genome during an im- 
mune response. 

Discussion 

Upon immune activation, DNA methylation in B cells is repro- 
grammed in a dramatic fashion at the GC stage and is not restored 
on differentiation to a memory or PC fate (Fig. 2B). The resulting 
similarity in DNA methylation patterns between the two post-GC 
cell types is unexpected, since the two cell types have distinct 
transcriptional profiles (Fig. 3). These results contrast to other 
differentiation models, where divergent cell fates derived from 



a multipotent/pluripotent progenitor adopt different epigenetic 
patterns (Mikkelsen et al. 2007; Ji et al. 2010). It is possible that 
a common methylation signature between memory and PCs par- 
tially explains the more rapid PC differentiation from memory 
B cells upon subsequent rechallenge with the same antigen. Our 
observations in the B-cell compartment indicate that both B and 
T cells (Messi et al. 2003; Kersh et al. 2006) employ epigenetic 
mechanisms in generating effective memory responses. 

While it is thought that repetitive elements are constitutively 
methylated in somatic cells, our results unexpectedly suggest that 
Alu elements display differential methylation patterns during im- 
mune activation within individuals (Figs. 3, 5). Our analysis also 
indicated enrichment of several other types of repeats by MIRA; 
however, we could not determine whether they are differentially 
methylated in B-cell subsets due to poor mapability at unique ge- 
nomic locations (data not shown). A substantial portion of the cell- 
type-specific differential methylation of Alus documented here 
correlates with decreased expression of the de novo methyl- 
transferase DNMT3A at the GC stage (Fig. 1C,D), consistent with 
a model delineated by Jones and colleagues Qeong et al. 2009). 
Diminished expression of DNMT3A, coupled with increased rep- 
licative demand could lead to replication-dependent loss of meth- 
ylation at repetitive elements. Further, the pattern of methylation 
loss at Alu elements is suggestive of an alteration in nucleosome 
positioning at the 5' and 3' ends of these elements, consistent with 
a global alteration in chromatin architecture. It is also conceivable 
that the DNA methylation changes observed here reflect a shift 
from homogenous to heterogeneous methylation patterning as 
previously described (Xie et al. 2011). 

The sheer number of significant changes in DNA methylation 
at Alus (14,000 individual^ luY elements) predicts that differential 
methylation in repetitive elements may contribute to differential 
gene expression within human populations in a manner analo- 
gous to that observed previously in animal models (Bultman et al. 
1994). In addition, it is likely that a subset of these methylation 
changes reflects or contributes to global changes in nuclear ar- 
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Figure 5. Widespread loss of DNA methylation at Alu elements in GC B cells. (A) Distribution of DMRs 
between naive and GC B cells at various genomic features (for definition of each genomic feature, see 
Methods). The bar graph illustrates the number of loss-of-methylation DMRs (naive > GC) in green and 
gain-of-methylation DMRs (naive < GC) in red that overlap with the indicated genomic feature. 
(B) Percentage of elements within each Alu subfamily (J, S, and Y) with loss (N > GC) or gain (N < GC) of 
methylation in GC compared to naive B cells by MIRA-seq analysis. (C) Validation of differentially 
methylated Alu elements in naive and GC B-cell populations by genomic bisulfite sequencing. The 
methylation status of CpGs across nine Alu elements was analyzed. After bisulfite treatment and PCR 
amplification of the Alus, at least eight individual clones of the PCR product were analyzed. The per- 
centage of methylated cytosine (/-axis) was calculated at individual CpG positions relative to the 5' end 
of each Alu element (x-axis). The genomic location (hgl 9) of the Alu elements analyzed are listed above 
each graph. 



DNA methylation upon immune acti- 
vation resembles the epigenetic repro- 
gramming events that occur during the 
generation of induced pluripotent stem 
(iPS) cells from somatic cells (Wernig 
et al. 2007; Mikkelsen et al. 2008; Kim 
et al.2010). DNA demethylation during 
this reprogramming process is a pre- 
requisite to achieve pluripotency status; 
inefficient DNA demethylation results 
in a partially de-differentiated state 
(Mikkelsen et al. 2008). Takahashi and 
Yamanaka (2006) elegantly demonstrated 
that the spectrum of transcriptional and 
epigenetic reprogramming required to 
generate pluripotent cells from adult 
precursors could be induced by specific 
TFs. We propose that the loss of DNA 
methylation we have described during 
the naive to GC B-cell transition permits 
GC B cells to acquire the potential to 
differentiate toward memory or PC fates, 
and to generate the differential response 
to antigenic challenge inherent to mem- 
ory cell populations. 

Methods 

Isolation of tonsillar B cells 

FACS purification of the B-cell subset 
was performed as described in Lai et al. 
(2010). Relevant clinical data of tonsil 
samples are included in Supplemental 
Table S8. The following antibodies were 
used: PE-IgD (BD Pharmingen), APC- 
CD38 (eBioscience, HIT2), and FITC- 
CD20 (eBioscience, 2H7). Propidium 
idodide was use to exclude dead cells. 
Sorting was performed using the BD FACS 
Aria II at the NIEHS Flow Cytometry 
Center. 



chitecture associated with the acquisition of a multipotent phe- 
notype. It will be essential to determine whether the methylation 
dynamics at Alu elements result in changes to their transcriptional 
status (Hellmann-Blumberg et al. 1993), as Alu transcripts have the 
capacity to modulate the global gene expression program through 
multiple mechanisms (Ponicsan et al. 2010). 

The GC B cell is a key intermediate in T-cell-dependent hu- 
moral immune responses. Like stem or progenitor cell populations, 
it is self-renewing, requires a specialized biological niche and sig- 
naling environment, and ultimately results in the formation of 
multiple, distinct cell types (Weissman 2000). Widespread loss of 



Nucleic acid extraction 
and manipulation 

Genomic DNA and total RNA were co- 
purified from purified cell types using the 
DNA/RNA Isolation Kit (Qiagen). RT-PCR 
was performed as described in Lai et al. 
(2010). For genomic bisulfite sequencing, 
the EZ DNA Methylation Gold Kit (Zymos) 
was used for bisulfite conversion of genomic DNA isolated from 
B-cell subsets. After PCR amplification, PCR products were cloned 
onto TOPO-XL vectors (Invitrogen) for subsequent sequencing. 



Immunohistochemistry and immunofluorescence 

Immunohistochemical and immunofluorescence staining of hu- 
man tonsil sections were performed as described in Fujita et al. 
(2004) and Jaye et al. (2007). Primary antibodies used are as fol- 
lows: DNMT1 (Novus Biologicals, NB100-264), DNMT3A (Santa 
Cruz Biotechnology, H295), BCL6 (Dako, PG-B6p), DNMT3B 
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Figure 6. Distribution of differentially methylated Alus across the genome. (A) Distance of differen- 
tially methylated AluYs to the nearest TSS. The distribution of all AluY (dashed purple line) and both types 
of differentially methylated Alus, loss-of-methylation (N > GC) and gain-of-methylation (N < GC), are 
shown. (B) Circos plot (Krzywinski et al. 2009) illustrating the density of differentially methylated AluY 
elements (middle ring, GC > naive; inner ring, naive > GC) on each chromosome. Color density (from 
yellow to red) represents the percentage of AluY elements (from 0% to 20%) in each 5-Mb bin that is 
differentially methylated (in the given direction). (C) Average normalized depth of mapped reads from 
MIRA-seq across a consensus profile for AluY; the position in base pair is shown on the x-axis. Percentage 
of AluYs with CpG at each position is also shown. (Left) The average depth of AluYs that have lost 
methylation in GC B cells; (right) the average depth of AluYs that have gained methylation in GC B cells 
compared to naive B cells. (D) Analysis of methylation status of genomic regions neighboring differ- 
entially methylated Alus by genomic bisulfite sequencing. Regions proximal to four differentially 
methylated Alus in Figure 5C were analyzed. After bisulfite treatment and PCR amplification, at least 
eight individual clones of the PCR product were analyzed. The percentage of methylated cytosine 
(/-axis) was calculated at individual CpG positions relative to the 5' end of each Alu element (x-axis). The 
genomic locations (hgl 9) of the Alu elements analyzed are listed above each graph. Methylation status 
of CpGs within individual PCR products are shown for five Alu elements in Supplemental Figure S4. 
(£) Overall DNA methylation dynamics at chromosomal regions with varying degrees of gene density 
between naive and GC B cells. We analyzed 5759 genomic bins of 0.5M bp for GC/naive adjusted RPKM 
ratios (there are a total of 61 46 bins in the genome, but only those bins with GC content between 30% 
and 60% were analyzed). Gene density for a bin is defined as the fraction of positions that are within at 
least one RefSeq gene body. The data were partitioned into regions with low gene density (gene density 
< 5%, N = 1 1 37, red line), intermediate gene density (5% < gene density < 40%, N = 1 533, green line), 
and high gene density (gene density > 40%, N = 31 00, purple line). 



(IMEGENEX, IMG-184A), MTA3 (Fujita 
et al. 2004), 5mC (Eurogentec, 33D3), 
and 5hmC (Active Motif, 59.1) 

Gene expression microarray 

Fifty nanograms of total RNA from B-cell 
subsets was used for cDNA synthesis and 
amplification using the WTA2 kit (Sigma). 
Amplified cDNA was labeled with Cy3 ran- 
dom nonomers (TriLink BioTechnologies) 
and hybridized onto NimbleGen Gene Ex- 
pression 12xl35K microarrays following 
the manufacturer's protocol (454 Life Sci- 
ences [Roche]). The microarray slides were 
scanned using the Agilent G2565BA 
DNA Microarray Scanner. Images were 
processed using the NimbleScan software. 
A three-step robust multi-array average 
(RMA) normalization (Ding et al. 2008) 
and batch correction (Johnson and Li 
2007) were performed on the data set. The 
RefSeq gene level expression data sets were 
generated by averaging expression signal 
from probes that overlap transcription 
start and end sites. Differentially expressed 
genes were defined using the criteria 
ANOVA P < 0.01 and fold change > 3. 

Global DNA methylation analysis 

MIRA-chip, data normalization, and iden- 
tification of methylation peaks were per- 
formed as described in Lai et al. (2010). 
NimbleGen 2.1M Deluxe Human Pro- 
moter Arrays (454 Life Sciences [Roche]) 
were used. Microarray probes are between 
50 and 75 bp in length. Detailed methods 
for peak calling, identification of DMRs, 
analysis of average methylation signals at 
genomic features, GSEA, and TF motif 
analysis are included in Supplemental 
Methods. For MIRA-seq, 1 fxg of genomic 
DNA from purified naive or GC B cells was 
used to perform MIRA for library prepara- 
tion. End repair, A-tailing, and adapter 
ligation steps were performed following 
Illumina's protocol. Adapter ligated DNA 
was size-selected for 250-600 bp on a 3% 
NuSieve 3:1 Agarose gel. Gel purified DNA 
(Qiagen) was then treated with sodium 
bisulfite for 5 h (EZ DNA Methylation- 
Gold Kit, Zymo Research), followed by 18 
cycles of PCR amplification using Pfu 
Turbo Cx (Stratagene). The libraries were 
sequenced on a Genome Analyzer IIx 
(Illumina) as paired-end 76mers; three 
lanes were run for each sample. Details of 
bioinformatics analysis at Alu elements 
are described in Supplemental Methods. 

Data access 

Microarray expression and DNA methyl- 
ation data have been submitted to the 
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NCBI Gene Expression Omnibus (GEO; http://www.ncbi.nlm. 
nih.gov/geo/) under accession number GSE24919. The MIRA-seq 
data have been submitted to the NCBI GEO under accession 
number GSE42386. 
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